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ABSTRACT 


Context. Supermassive stars and quasi-stars (massive stars with a central black hole) are both considered as potential progenitors for 
the formation of supermassive black holes. They are expected to form from rapidly accreting protostars in massive primordial halos. 

Aims. We explore how long rapidly accreting protostars remain on the Hayashi track, implying large protostellar radii and weak 
accretion luminosity feedback. We assess the potential role of energy production in the nuclear core, and determine what regulates 
the evolution of such protostars into quasi-stars or supermassive stars. 

Methods. We followed the contraction of characteristic mass shells in rapidly accreting protostars, and inferred the timescales for 
them to reach nuclear densities. We compared the characteristic timescales for nuclear burning with those for which the extended 
protostellar envelope can be maintained. 

Results. We find that the extended envelope can be maintained up to protostellar masses of 3.6 x 10871? Mo, where rn denotes 
the accretion rate in solar masses per year. We expect the nuclear core to exhaust its hydrogen content in 7 x 10° yr. If accretion 
rates m > 0.14 can still be maintained at this point, a black hole may form within the accreting envelope, leading to a quasi-star. 
Alternatively, the accreting object will gravitationally contract to become a main-sequence supermassive star. 

Conclusions. Due to the limited gas reservoir in typical 10” Mo dark matter halos, the accretion rate onto the central object may drop 
at late times, implying the formation of supermassive stars as the typical outcome of direct collapse. However, if high accretion rates 
are maintained, a quasi-star with an interior black hole may form. 


Key words. cosmology: theory — dark ages, reionization, first stars — stars: protostars — stars: formation 


1. Introduction 


Supermassive black holes (SMBHs) with more than 10° Mo 
have been observed at z > 6 (Fan et al. 2001, 2003, 2004, 2006), 
and recently even at z = 7.085 (Mortlock et al. 2011). While the- 
oretical considerations indicate a variety of potential pathways, 
including direct collapse to a black hole, collapse to a supermas- 
sive star, collapse of stellar clusters or clusters of black holes 
(Rees 1984; Haiman 2006; Volonteri & Bellovary 2012), sub- 
stantial efforts have been made to assess the validity of such 
scenarios. 


1.1. Black hole formation scenarios 


As the most straightforward possibility, one could imagine these 
SMBHs to originate from remnants of the first stars, which were 
proposed to be very massive with 100—1000 Mo (Abel et al. 
2002; Bromm & Larson 2004). Stellar evolution calculations by 
Heger & Woosley (2002) indeed show that black holes may form 
for stellar masses between 10-30 Mo, or 100—300 Mo, thus 
yielding the seeds for further growth to SMBHs. 


Article published by EDP Sciences 


New calculations following the evolution of primordial 
star formation beyond the first core, however, show that 
fragmentation can occur efficiently (Stacy et al. 2010, 2012; 
Turk et al. 2009; Clark et al. 2011; Greif et al. 2011; Smith 
et al. 2011). While the effect of the accretion luminosity on the 
protostellar core seems to be minor (Smith et al. 2011, 2012), 
UV feedback during the main sequence phase appears to set 
the upper limit to the stellar mass to 50—100 Mo (Hosokawa 
et al. 2011, 2012b; Susa 2013). Even if very massive stars could 
form, the resulting HII region is likely to inhibit accretion for at 
least 108 yrs (Milosavljević et al. 2009a,b). An additional prob- 
lem has been identified by Whalen & Fryer (2012), who showed 
that black holes resulting from such stars would be born with 
strong dynamical kicks, implying that they should readily be 
expelled from the halo. 


On the other hand, dense stellar clusters may form at high 
redshift in the presence of trace amounts of dust (Schneider 
et al. 2003, 2008, 2012a,b; Omukai et al. 2005, 2012; Clark 
et al. 2008; Dopcke et al. 2011, 2013; Klessen et al. 2012). The 
relativistic instability of such clusters has been explored in a 
semi-analytic framework by Devecchi & Volonteri (2009) and 
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Devecchi et al. (2010, 2012), finding characteristic black hole 
masses of up to ~3000 Mo. The latter is already substantially 
higher than the masses of the first stars, though it may still be dif- 
ficult to obtain masses as high as 10? Mo at z = 7. In the presence 
of a significant black hole spin, which will be rapidly obtained in 
the presence of efficient accretion, one typically expects an addi- 
tional growth by about four orders of magnitude, thus requiring 
seed masses on the order of 10° Mo (Shapiro 2005). 

The formation of such massive seeds has been consid- 
ered in the context of the direct collapse model, where gas in 
a 107—108 Mg halo is expected to collapse without fragmentation 
into a single central object (Koushiappas et al. 2004; Begelman 
et al. 2006; Lodato & Natarajan 2006; Spaans & Silk 2006; 
Volonteri et al. 2008; Begelman & Shlosman 2009). In order to 
avoid fragmentation, molecular hydrogen needs to be efficiently 
dissociated, requiring a strong ambient UV field (Omukai 2001; 
Bromm & Loeb 2003; Shang et al. 2010; Schleicher et al. 2010; 
Latif et al. 2011; Van Borm & Spaans 2013). Such radiation 
backgrounds may have been provided by nearby starburst galax- 
ies, which may indeed occur frequently enough to explain the 
observed abundance of SMBHs at z ~ 6 (Dijkstra et al. 2008; 
Agarwal et al. 2012). As an alternative, Inayoshi & Omukai 
(2012) suggested that black hole formation may occur in the 
presence of strong shocks, leading to collisional dissociation of 
molecular hydrogen. 

Using numerical simulations, Wise et al. (2008) have mod- 
eled the gravitational collapse of massive primordial halos cool- 
ing via atomic hydrogen, reporting an isothermal density profile 
and angular momentum transport via bar-like instabilities dur- 
ing the formation of the first core. Regan & Haehnelt (2009) 
presented the first study following the evolution beyond the first 
core, and modeling the formation of self-gravitating disks on 
parsec-scales. While these studies employed a typical resolu- 
tion of 16 cells per Jeans length, it was recently demonstrated 
that turbulent structures can only be resolved with a resolution 
of at least 32 cells per Jeans length, preferably more (Sur et al. 
2010; Federrath et al. 2011; Turk et al. 2012). In a set of high- 
resolution simulations following the formation of the first core, 
Latif et al. (2013a,b) found extended turbulent structures in the 
center of massive primordial halos, but no signs of simple disks 
or bar-like instabilities. Disks do, however, form at later stages of 
the evolution, with masses of ~1000 Mo on scales of 30 AU, and 
characteristic accretion rates of ~1 Mo yr™! (Latif et al. 2013c). 
As a result, very massive central objects may indeed form in 
rather short cosmic times. The dynamics at early times have also 
been explored by Choi et al. (2013), finding a somewhat differ- 
ent result including the formation of toroidal structures on scales 
of several parsec. It is, however, not fully clear to which extent 
the latter is a result of the simplified initial conditions employed 
in their calculation. To provide more quantitative predictions re- 
garding the conditions where black holes may form, Prieto et al. 
(2013) explored the correlations of the baryon spin with the spin 
of the dark matter, showing that the resulting correlation is, how- 
ever, weak, and that the baryon properties cannot be naively ex- 
trapolated from the dark matter. In a similar study, de Souza 
et al. (2013) report that the dynamical properties of the dark 
matter halo may leave an imprint on the resulting stellar mass. 


1.2. Quasi-stars vs. supermassive stars 


While the accumulation of high masses seems feasible from a 
hydrodynamical point of view, the resulting object is still consid- 
erably more uncertain. As a first step, we expect the formation 
of a massive protostar, as the gas becomes optically thick at high 
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densities (Omukai & Palla 2001, 2003; Hosokawa et al. 2012a). 
It is, however, unclear how these objects are going to evolve, and 
whether they will form a supermassive star or a quasi-star. 


In the present context, we define a supermassive main se- 
quence star as a conventional star with very high mass (M = 
10°—10° Mo) (see also Shapiro & Teukolsky 1986), which pre- 
dominantly generates its energy via nuclear burning. A quasi- 
star refers to an object with mass similar to that of a super- 
massive star, but whose central core has collapsed into a black 
hole (Begelman et al. 2006; Begelman 2010). The progenitors of 
both types of objects are rapidly accreting high-mass protostars, 
and we aim to explore in this paper under which conditions they 
evolve into a quasi-star or a supermassive main-sequence star. 


To explore the evolution of such protostars, Hosokawa et al. 
(2012a) recently followed the evolution of rapidly accreting pro- 
tostars, showing that they expand as cool supergiants, thus in- 
hibiting the feedback from accretion luminosity and potentially 
allowing accretion to proceed for very long times. Hosokawa 
et al. also found that nuclear burning starts at stellar masses of 
about ~50 Mo, and is maintained until the end of their calcu- 
lation at 10° Mo. Therefore, their model indicates no transition 
toward a quasi-star at least during the evolutionary phase con- 
sidered. It is, however, important to assess how long this phase 
of efficient accretion can be maintained, and under which con- 
ditions a transition to a supermassive main-sequence star can 
be expected, which can potentially influence the accretion flow 
via radiative feedback (Omukai & Inutsuka 2002; Johnson et al. 
2011, 2012). While the majority of such supermassive stars di- 
rectly collapses into a black hole (Fryer & Heger 2011), a small 
mass window exists around 55 000 Mo where violent supernova 
explosions of up to 10% erg may occur (Johnson et al. 2013). 


Following a different approach, Begelman et al. (2006) and 
Begelman (2010) have proposed the formation of a black hole in 
the interior of rapidly accreting objects, leading to quasi-stars as 
the progenitors of SMBHs. Begelman et al. (2006) argued that 
for quasi-stars with 10*-10° Mo, the typical accretion timescale 
is considerably shorter than the timescale for nuclear burning, 
and as a result, the latter may not be able to stop the col- 
lapse of the central core. A potential advantage of such a con- 
figuration is that the accretion onto the central object is not 
limited by the Eddington accretion rate of the black hole, but 
by the Eddington accretion rate of the more massive quasi- 
star (Begelman et al. 2008). Begelman (2010) considers objects 
with more than 10° Mo, which first evolve toward the main se- 
quence as a supermassive star, but then form a black hole af- 
ter an extended phase of nuclear burning. Ball et al. (2011) fol- 
lowed the evolution of black holes in a quasi-star employing the 
Cambridge STARS stellar evolution package (Eggleton 1971; 
Pols et al. 1995), finding that the black hole may efficiently 
accrete 10% of the stellar mass before hydrostatic equilibrium 
breaks down. They also report that the results are sensitive 
to the boundary condition in the interior. The contraction of 
an embedded isothermal core in a stellar envelope was fur- 
ther explored by Ball et al. (2012) in the framework of the 
Schénberg-Chandrasekhar limit. 

While the models for the quasi-stars typically employed a 
Thomson-scattering opacity, Hosokawa et al. (2012a) reported 
that H` rather than Thomson scattering dominates the opacity in 
the protostellar atmosphere. They find luminosities close to the 
Eddington luminosity, 


4rGMmpc 


OT 


M 
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where G is the gravitational constant, mp the proton mass, c the 
speed of light, or the Thomson scattering cross section and M 
the mass of the star. The temperatures in the atmosphere are 
rather cool, ~5000 K, as the protostars are on the Hayashi track. 
As a result, the characteristic radii evolve as 


1/2 
M 
Rin = 2.6 X 10° lé . (2) 
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The protostars are thus considerably more extended and show an 
increasing stellar radius as a function of mass, while the models 
of Begelman et al. (2006) indicate a constant radius as a function 
of stellar mass. Begelman (2010) consider a mass-dependent 
stellar radius, which is, however, still more compact due to the 
assumption of Thomson scattering opacities. 

The difference between these scenarios may, however, turn 
out to be highly relevant, as the stellar radius regulates the tem- 
perature on the surface, and thus the strength of protostellar feed- 
back. In fact, the behavior reported by Hosokawa et al. (2012a) 
is well-known (Omukai & Palla 2001, 2003) in cases where the 
accretion timescale, 


M 
lacc = =; 3 
P (3) 


is much smaller than the Kelvin-Helmholtz timescale, 


(4) 


However, Hosokawa et al. (2012a) report that this behav- 
ior extends into the regime where Ten < facc, i.e. where 
Kelvin-Helmholtz contraction is faster than mass growth via ac- 
cretion. The physical mechanism which allows this phase to con- 
tinue will be discussed in detail in Sect. 2. Taking the maximum 
luminosity of the star 


11/2 -1/2 
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which assumes Kramers’ opacity x œ of 27 (Hayashi et al. 
1962) and an initial radius as given in Eq. (2), they show that 
both timescales are equal at protostellar masses of 
Re" 6 
or (om) Ge 


where we parametrized the accretion rate as M = m Mg yr™!. On 
the other hand, the accreting star remains on the Hayashi track 
up to stellar masses of 1000 Mo, the highest mass reached in 
the calculation of Hosokawa et al. (2012a). The extended enve- 
lope locked at Teg ~ 5000 K allows efficient accretion with only 
moderate feedback for a longer period. 

In this paper, we aim to assess how long this efficient ac- 
cretion phase can be maintained beyond M ~ 1000 Mo without 
strong feedback from the protostar. In addition, we will discuss 
the potential impact of nuclear burning for the evolution of the 
accreting objects, and the conditions under which a quasi-star as 
opposed to a supermassive star may form. For this purpose, we 
explore the interplay of mass accretion and Kelvin-Helmholtz 
contraction in Sect. 2. We calculate the impact of nuclear burn- 
ing in Sect. 3. A final discussion of our results is provided in 
Sect. 4. The impact of additional processes such as deuterium 
shell burning and hydrogen burning via the pp-chain is assessed 
in the Appendix. 


2. The interplay of mass accretion 
and Kelvin-Helmholtz contraction 


In the following, we will sketch why the protostar may main- 
tain large envelopes even if the Kelvin-Helmholtz timescale be- 
comes shorter than the timescale for accretion. For this pur- 
pose, we consider mass shells of enclosed mass M, located at 
radii R(M, f). A mass shell M forms at the time M/M, with an 
initial radius as given by Eq. (2). The radii of these mass shells 
evolve on their Kelvin-Helmholtz timescale, given as 


tku(M, R) = EE (7) 
R(M, t) Lead 

As a result, we have 

dR RMD me 

dr tku(M, R) 


The equation is integrated between the initial radius Rin when 
the mass shell forms, given in (2), corresponding to the 
time Gol = M/M, to the current radius R at time t. The 
integration yields 

1 1 4am 

= = — + 

R Rn o1rM 


(t — tin(M))- (9) 


Inserting the initial radius as given in Eq. (2) and expressing the 
result in astrophysical units, we obtain 


1000 Ro ` 1000 Ro 
R—-2.6 x 103(M/Mo)!/2 
1.16 x 107! yr! (t — tin(M 
S 6x10 yr (t ( 1 (10) 
M/Mo 


The first term on the right-hand side thus dominates right af- 
ter the formation of a mass shell and determines its initial 
radius, while the second term describes its evolution due to 
Kelvin-Helmholtz contraction. It is remarkable that at late times, 
the second term dominates and the evolution of the radii ap- 
pears to be independent of their initial position. This can be un- 
derstood, as our expression for the Kelvin-Helmholtz timescale, 
Eq. (7), scales with the inverse radius of the mass shell. As a re- 
sult, the evolution slows down during the contraction, implying 
that the Kelvin-Helmholtz timescale may increase significantly 
in the interior. An example for the resulting structure in a 10° Mo 
star is given in Fig. 1, including both the regime where the first 
and the second term of the equation dominate. 

We now aim to estimate when a given mass shell reaches the 
densities of nuclear burning, which we take as pmc ~ 1 g cm~? 
(Hosokawa et al. 2012a). We obtain a radius of nuclear burning, 
which is given as 


1/3 1/3 
M 

) S L2ro( 2) . 
M, 
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R 3M 
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We note that the above provides only a rough estimate of the 
radius of the nuclear core. If the latter were more compact, the 
timescale for a given mass scale to reach nuclear burning would 
become even longer, implying that more mass could be accreted 
during that time. A comparison with Eq. (2) shows that this ra- 
dius always remains smaller than the radius of the star by at least 
a factor of 1000, and in fact increases more gradually with stel- 
lar mass. A given mass shell M will only reach nuclear densi- 
ties on timescales much longer than the initial Kelvin-Helmholtz 
timescale tgu, which increases as R~!. We can thus determine the 


(11) 
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Fig. 1. Enclosed mass as a function of radius in the protostar calcu- 
lated from Eq. (10), assuming a total stellar mass of 10° Mọ and 
mass accretion rates m = 0.1,1.,10. The current time rie calculated 
as t = (M/M,)/m. For a given mass shell, higher accretion rates imply 
larger radii, as the star had less time to contract. 


time when nuclear densities are reached by equating (11) with 
the second term in (10). As a result, we obtain 


Vun 
At = t — tin(M) = 710 cl) . (12) 
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During the time interval Ar, the star will accrete an additional 


mass AM = MAt. Considering a given mass shell, we can 
calculate the ratio 

AM _ MAt m\'? 

— = — ~710m|—]_, (13) 
M M Mo 


which describes how much additional mass is accreted before 
the shell reaches nuclear densities. The ratio drops below 1 only 
for protostellar masses of 


M > 3.6 x 108 m? Mo. (14) 


While for typical accretion rates of rn ~ 107°, this happens al- 
ready at mass scales of order unity, this transition occurs only at 
substantially higher masses above 1000 Mo for m > 1077. While 
a given shell evolves toward nuclear densities, the amount of ad- 
ditional mass that is accreted is thus substantially higher than the 
mass in that shell. From Eq. (12), we can further derive the mass 
in the nuclear core in the limit that t > tin( M). We obtain 


, ii 

Myuc = Lil Mo = oe Mo, 
where we defined tun = t/(710 yr). As we neglected the 
term M/M in Eq. (12), we note that the derivation here implic- 
itly assumes a sufficient mass supply to the protostar in order to 
feed the core. Clearly, this approximation will break down at the 
mass scale derived in (14). At that point, the protostellar mass 
may become approximately constant, and we expect a transition 
toward supermassive main sequence stars. 

From this model, we already obtain a set of relevant con- 
clusions concerning the evolution of the protostar. In particular, 
if the accretion rate is high, a substantial amount of additional 
matter is accreted before a given shell reaches nuclear densities. 
As a result, the evolution of the protostar is dominated not by 
the interior shells, but by the additional matter which is accreted 


(15) 
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during the contraction. This behavior will only change when 
the timescale of accretion becomes comparable to the timescale 
on which the outer mass shell is able to reach nuclear densi- 
ties. As we show here, the latter is significantly larger than the 
Kelvin-Helmholtz timescale of the star, and scales as R,!.. The 
transition occurs when a critical mass scale of 3.6 x 10°? Mo 
is reached, implying a transition toward a supermassive main se- 
quence star. We note that the transition may even occur earlier 
if m decreases with time. 


3. The nuclear evolution 


In this section, we will assess the potential role of nuclear 
burning. For this purpose, we show that the transition to the 
CNO-cycle rapidly occurs and regulates the production of he- 
lium in the nuclear core. We will then calculate under which 
conditions the fuel in the nuclear core will be exhausted while 
the protostar maintains its extended envelope. The latter implies 
the potential formation of a black hole in the interior, thus a 
transition into a quasi-star. 


3.1. Importance of the CNO-cycle 


While the pp-cycle will dominate in the very beginning (see 
appendix), we show here that the transition to the CNO-cycle 
rapidly occurs as a result of helium burning via the triple-a 
process. Following Padmanabhan (2000), the energy production 
rate via triple-a@ is given as 


5.1 x 1085Y? 


7 e744027/T9 erg ST el, (16) 


Bq 


where Y ~ 0.25 denotes the mass fraction of helium and To de- 
notes the gas temperature in units of 10° K. It is straightforward 
to check that the energy production rate changes by several or- 
ders of magnitude for temperature changes on the order of 10%. 
As shown by Hosokawa et al. (2012a), the highest temperatures 
with Tọ ~ 0.15 are expected in the center of the nuclear core, 
which has an enhanced nuclear density of 10 g cm™*. With a 
volume filling factor ej, ~ 0.01—0.1, the contribution of that 
region to helium burning is significant due to the steep temper- 
ature dependence. For the following estimates, we will assume 
that efficient mixing occurs throughout the nuclear core, imply- 
ing that the heavy elements produced in the central region will 
be available throughout the core. On the other hand, if mixing 
is inefficient, nuclear burning via the CNO cycle will dominate 
even more in the central region due to its higher metallicity. 

Noting that the energy released by a single triple-@ reaction 
corresponds to 1.166 x 1075 erg, the number of heavy elements 
produced per unit time is given as 


E3a Efill Myuc 


N3q = 1.166 x 10-5 erg. (17) 
The CNO mass in the core thus evolves as 

Mcno = (inte X Noo (18) 
and an integration yields 

Meno = 7.1 10°" ein Ga Mo. (19) 


Assuming efficient mixing, the resulting metallicity in the 
nuclear core is then 


Mcno 


Z= = 6.3 X 10°? t710 emm, (20) 


nuc 
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We recall that the energy production rate in the CNO cycle is 
given as (Padmanabhan 2000) 


2/3 
T 


ECNO = erg g`! sl. (21) 


Since the energy production rate depends sensitively on the tem- 
perature, we expect its contribution in the innermost core to be 
dominant. For the temperature Ty = 0.15, it is straightforward to 
show that energy production via the CNO cycle becomes com- 
parable to the pp cycle for a metallicity of Ze = 2 x 10°. 
A comparison with Eq. (20) shows that the CNO cycle thus 
dominates after a short time of 


tono = Zen A6, (22) 


where we introduced ej = O.lég,-1. We note here that this 
timescale is probably not accurate, as our assumptions in Sect. 2 
(in particular concerning the Eddington luminosity) only be- 
come valid at later times. Nevertheless, this result implies that 
the CNO cycle can be expected to be relevant early on, and 
thus needs to be considered. As the production of elements be- 
comes increasingly efficient at late times, we expect that the ex- 
pression (20) will yield a reasonable estimate in the period of 
interest. 

With the above assumptions, the CNO luminosity due to 
nuclear burning is given as 


Z 
Leno = EcNo.Z. E Maue en, (23) 


Ze 
where €cno,z, denotes the energy production rate via the 
CNO cycle evaluated at the critical metallicity Z.. Inserting (20), 
we obtain 


Leno = 1.3 x 10° 2, Eis Lo. 


u 510 (24) 


The Eddington luminosity of the nuclear core is given as 


Ledd,core = 3.8 x ri" is 


710 Lo. 


(25) 


The ratio of these luminosities to the Eddington luminosity of 
the star is given in Fig. 2. Even here, the luminosity produced 
by nuclear burning approaches the Eddington luminosity of the 
star only around stellar masses of 108 Mo, implying no rele- 
vant impact on the stellar evolution during the earlier stages. A 
comparison of these expressions yields the timescale 


te = 2.1 x 10'e 4 YE (26) 


For enn 1 ~ 1, as indicated by Hosokawa et al. (2012a) for m = 
0.1, the CNO luminosity would exceed the Eddington luminosity 
of the core after ~2 x 10* yr, shortly after the end of their sim- 
ulations. However, when the Eddington luminosity is reached, 
one would expect an expansion of the nuclear core, implying 
lower densities and an adiabatically decreasing temperature. Due 
to this thermostat, it is likely that the core will adjust to a state 
maintaining its Eddington luminosity. 

We will now demonstrate that the nuclear core is not con- 
verted to helium before the critical timescale te after which 
the core radiates at its Eddington luminosity. We consider the 
production of helium by the CNO process, given as 


Leno 


Keeseren 27 
4.3 x 10° erg (27) 


Nu cu = 


and the helium mass production rate 


Mye,cno = 12mp X Nos co, (28) 
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Fig. 2. Ratio of nuclear luminosity provided via the CNO-cycle over the 
Eddington luminosity of the star as a function of protostellar mass for 
different mass accretion rates. The CNO luminosity is calculated here as 
the minimum of the estimate in Eq. (24) and the Eddington luminosity 
of the nuclear core in Eq. (25). At high protostellar masses, one can 
clearly recognize the transition to the regime where nuclear burning is 
limited by the Eddington luminosity of the core. Higher accretion rates 
imply that less time was available for contraction, thus reducing the 
nuclear luminosities. 


An integration yields 


7/2 


Mue,cno = 1.0 x II een Gi Mo. (29) 


A comparison with the mass of the core indicates that they 
become comparable at 


te,cNno = 2.25 x 10°54) Y5 (30) 


implying that a pure helium core will only be formed once the 
nuclear luminosity is equal to the Eddington luminosity of the 
core. 


3.2. Formation of a helium core 


We will now calculate the helium production rate under the as- 
sumption that the formation of a helium core occurs once the 
nuclear luminosity equals the Eddington luminosity of the core. 
Thus, 


. LEda,core 
Nie = E (31) 
The helium mass production rate is thus given as 

Mue = 4m, X Nue (32) 
and an integration yields 

Mne = 1.05 x 1074 ia Mo. (33) 


Equating Eq. (33) with Eq. (15), we find that the nuclear core 
becomes a pure helium core after ty ~ 6.8 x 10° yr. We note 
that at this point, the mass of the nuclear core is given as My = 
9.4 x 10° Mo. In order to illustrate how the latter translates 
into total stellar masses, the ratio of helium mass to total mass 
in the core is displayed in Fig. 3 for different accretion rates, 
showing that the hydrogen fuel is indeed exhausted for stel- 
lar masses of 3 x 10°-3 x 10’ Mo depending on the accretion 
rate. After that point, the burning of heavier elements will oc- 
cur, but only last for a short period. As the mass of the core is 
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Fig. 3. Fraction of helium in the nuclear core, corresponding to the 
amount of exhausted fuel for nuclear burning, as a function of proto- 
stellar mass for different accretion rates. Higher accretion rates imply 
less time for protostellar contraction. As a result, nuclear densities are 
reached later, implying less time for nuclear burning. 


considerably higher than the Chandrasekhar mass scale and the 
Tolman-Oppenheimer- Volkoff mass scale, the core may collapse 
and form a black hole within the star, a so-called quasi-star, as 
proposed by Begelman et al. (2006) and Begelman (2010). The 
interior black hole may then continue accreting from the stellar 
envelope, while the accretion rate of the quasi-star will only be 
limited by the Eddington rate corresponding to the total mass of 
the configuration. The derivation of t and My, however, implic- 
itly assumes that enough mass has been accreted in order to build 
up the core, requiring that 
My < tc. (34) 
To form a quasi-star, we thus require an accretion rate of M > 
M f /t fo Le. 
m > 0.14. (35) 
For lower accretion rates, we expect that the extended envelope 
of the star cannot be maintained at late times, implying that the 
star will contract and form a supermassive main sequence star. 
During this period, both the stellar mass and the mass in the nu- 
clear core will be approximately constant, implying a character- 
istic hydrogen burning timescale of 2 x 10° yr. For stellar masses 


above 300 Mo, these are expected to collapse to massive black 
holes (Fryer et al. 2001; Fryer & Heger 2011). 


4. Discussion and outlook 


To investigate the interplay of accretion and Kelvin-Helmholtz- 
contraction of the protostar, we have presented an analytical 
model following the evolution of the mass shells from the ex- 
tended envelope down to nuclear densities. Based on this model, 
we are able to explain the fundamental result by Hosokawa 
et al. (2012a) that the accreting protostar keeps its extended 
envelope beyond the adiabatic accretion phase. The reason is 
that the timescale on which a given mass shell may reach 
nuclear densities is much longer than the Kelvin-Helmholtz 
timescale of the star, as the protostellar radius is consider- 
ably larger than the radius of the nuclear core. Instead of the 
Kelvin-Helmholtz timescale evaluated at the protostellar radius, 
the relevant timescale is thus the Kelvin-Helmholtz timescale 
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Table 1. Pathways to the formation of main-sequence supermassive 
stars and quasi-stars depending on the critical mass scale MA) = 
3.6 X 108 m? Mo and the accretion rate sn in solar masses per year. 


M < M.,(m) M > M.,(m), 
ms 0.14 
transition to 
main-sequence 


supermassive star 


M > Man 
m > 0.14 
collapse of the 
nuclear core, 
quasi-star 


protostar with 
extended envelope, 
cool atmospheres 


at the expected core radius after contraction, providing an ap- 
propriate estimate of the evolutionary timescale. As a conse- 
quence, rapidly accreting protostars are expected to keep their 
extended envelopes until stellar masses of 3.6 x 108m? Mo. 
At this point, the accretion timescale becomes longer than 
the Kelvin-Helmholtz timescale in the stellar interior, imply- 
ing that the star will subsequently contract toward the main se- 
quence. Then, we expect it to follow the typical evolution of a 
supermassive star as sketched by Shapiro & Teukolsky (1986). 


We further assessed the impact of nuclear burning for the 
evolution of these stars. While deuterium burning will never 
be relevant during the protostellar evolution, hydrogen burning 
via the pp-cycle could potentially affect the stellar evolution af- 
ter 10!5 years, which is, however, practically irrelevant. Due to 
the efficient metal production via the triple-@ process, a transi- 
tion to the CNO cycle is expected early on. The energy produc- 
tion by the CNO cycle is considerably more efficient, implying 
that the core can radiate at its Eddington luminosity. As a result, 
a helium core forms after 7 x 10° yr. The expected mass of the 
core is then 10° Mo, which provides an upper limit on the ini- 
tial mass of the interior black hole. Therefore, if the available 
gas reservoir is sufficient to maintain accretion rates m > 0.14 
until this point, the central core may collapse during the pro- 
tostellar evolution phase, and the object becomes a quasi-star. 
Alternatively, the protostar will evolve toward the main sequence 
to become a supermassive star. Such an object is then expected 
to collapse into a SMBH at the end of its lifetime (Fryer & 
Heger 2011). A summary of our results regarding the evolution 
of rapidly accreting protostars is given in Table 1. 


Our results have been derived based on a comparison of the 
mass accumulated by accretion with the characteristic mass scale 
where we expect a pure helium core. This mass scale is based 
on the reasonable assumption that the luminosity of the core is 
given by its Eddington luminosity, in agreement with the nu- 
merical results by Hosokawa et al. (2012a). Uncertainties are, 
however, present, as our calculation for the evolution of the core 
mass assumes that the luminosity of each mass shell is given by 
the Eddington luminosity. Following Hosokawa et al. (2012a), 
the latter is a good but not precise approximation. Further cor- 
rections can be expected in particular when approaching nuclear 
densities, when feedback from nuclear burning becomes signif- 
icant. We further note that in the regime where deviations from 
the hydrostatic equilibrium are significant, as during the collapse 
of the nuclear core or during the accretion onto the protostellar 
envelope, it is conceivable that the luminosity released from the 
gravitational potential may locally exceed the Eddington lumi- 
nosity if the accretion rate is high enough. Due to the momen- 
tum of the flow, we, however, expect that such accretion rates 
could still overcome the radiative feedback and effectively trap 
the photons in the optically thick material. We further note that 
Dotan & Shaviv (2012) propose a scenario where porosity in- 
duced by instabilities may lead to super-Eddington luminosities. 
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The results obtained differ from the scenario proposed by 
Begelman et al. (2006) and Begelman (2010), as we make use 
of new results on the protostellar evolution that were obtained 
in the recent study by Hosokawa et al. (2012a). In particu- 
lar, we note that Begelman et al. (2006) still assumed a time- 
independent radius of the quasi-stars, and both Begelman et al. 
(2006) and Begelman (2010) adopted Thomson scattering as the 
only opacity source. However, Hosokawa et al. (2012a) have 
shown that H" may significantly contribute to the opacity in the 
atmosphere. By considering the timescale required for a given 
mass shell to reach nuclear densities, we further have shown here 
that the star may remain on the Hayashi track for a much longer 
period than expected before. The latter gives rise to consider- 
able differences in the stellar evolution, and allows for significant 
amounts of accretion while feedback is weak. 


Additional effects could be introduced as a result of rotation, 
which was not included in this model. We thus expect our re- 
sults to provide an order of magnitude estimate for the critical 
accretion rate, which may be determined more accurately em- 
ploying stellar evolution calculations, but also a more realistic 
model for the time-dependent accretion rate. We also note here 
that Hosokawa et al. (2012a) reported a slight dependence of 
their results on the employed boundary conditions for the pro- 
tostar. While their main calculations adopted shock boundary 
conditions, they also explored the effect of lower-entropy ac- 
cretion provided by photospheric boundary conditions. In that 
case, a slightly higher accretion rate of 0.3 Mo yr! appears to 
be required to maintain the extended envelopes, but the overall 
evolution remains very similar. 


For comparison, we note that recent simulations by Latif 
et al. (2013c) reported accretion rates of ~1 Mo yr™! in halos 
with ~10’ Mo. Adopting a constant accretion rate, the transition 
to a supermassive main sequence star would occur at a stellar 
mass of 3.6 x 108 Mo. As the nuclear fuel is, however, exhausted 
earlier, a central black hole may form before this transition, giv- 
ing rise to an evolution as sketched by Ball et al. (2011). Due 
to the limited gas reservoir, and as here the halo mass is in fact 
comparable to the mass scale of the corresponding quasi-star, it 
seems likely that the accretion rate will substantially decrease at 
late times, implying that the transition toward the supermassive 
star should occur at an earlier stage which is then determined by 
the time evolution of the accretion rate, and that supermassive 
stars with 104—105 Mo might be the most generic outcome of 
the collapse. More massive objects appear to be possible at least 
in principle, if a larger gas reservoir is available in more massive 
dark matter halos. 


In order to assess the potential influence of rotation, 
we estimated the amount of rotational support in the cen- 
tral 10° Mọ clumps reported by Latif et al. (2013c), which 
varied between 5—20% in different simulations. With densities 
of ~107! g cm?, these are still orders of magnitudes below the 
characteristic densities within the protostar, such that we cannot 
yet draw strong conclusions regarding the final amount of rota- 
tional energy. However, we note that Stacy et al. (2012) reported 
a significant amount of rotation in primordial protostars based 
on numerical simulations of Greif et al. (2012), and a quite sim- 
ilar case can be expected here. As recently shown by Reisswig 
et al. (2013), the latter may have a substantial impact on the col- 
lapse of supermassive stars, implying the potential formation of 
a black hole binary and a subsequent merger, accompanied with 
efficient emission of gravitational waves. The latter provides a 
potential pathway of probing black hole formation scenarios 


with LISA!. In a narrow mass range around ~55 000 Mo, one 
may further expect the occurence of highly energetic supernovae 
with energies up to 10“ erg, which can be potentially detected 
with JWST? (Johnson et al. 2013). 

Returning to the fate of supermassive stars, their evolution 
in the presence of UV feedback was assessed by Omukai & 
Inutsuka (2002) and Johnson et al. (2012) in the case of spherical 
symmetry. These authors find that UV feedback is unable to stop 
accretion for rates above ~0.1 Mo yr™!. As a result, the expected 
stellar mass is 


m IT (36) 


3 

Muy ~ 10° Mo (5 
For accretion rates of ~0.1 Mo yrl, our results exceed this 
value, as the protostars remain on the Hayashi track up to a mass 
of ~3 x 10° Mo, making feedback inefficient. Therefore, taking 
protostellar evolution into account favors the formation of more 
massive objects. In all cases, a critical question concerns the time 
evolution of the accretion rate, since the transition to supermas- 
sive stars is regulated by the accretion rate at late times. While 
this paper provides a first assessment for the case of constant ac- 
cretion rates and spherical symmetry, the potential implications 
of time-dependent accretion rates need to be addressed in more 
detail in the future, along with the implications of rotation during 
protostellar evolution. 
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Appendix A: Additional nuclear processes 


In Sect. 3.1, we have shown that the transition to the CNO cycle 
rapidly occurs, implying that the latter will regulate the forma- 
tion of a helium core as discussed in Sect. 3.2. However, addi- 
tional nuclear processes are expected to simultaneously occur, 
which we discuss here for completeness. In the two sub-sections 
below, we assess the role of deuterium shell burning as well as 
the impact of hydrogen burning via the pp-chain during the stel- 
lar evolution, showing that these processes will only have a mi- 
nor impact on the stellar evolution. 


A.1. Deuterium burning 


As reported by Hosokawa & Omukai (2009), deuterium shell 
burning may occur even before nuclear burning starts in the core. 
The luminosity from deuterium shell burning is given as 


[D/H] l (A.D 


0.1/\2.5 x 10-5 


Lp = 15x 10° Gë 
where [D/H] denotes the deuterium abundance relative to hydro- 
gen. For an accretion rate of m = 0.1, Hosokawa et al. (2012a) 
demonstrated that deuterium shell burning starts only at proto- 
stellar masses of ~80 Mo, and presumably even higher masses in 
case of higher accretion rates. A comparison with the Eddington 
accretion rate in Eq. (1) thus shows that the luminosity from deu- 
terium shell burning will never become dominant in this mass 


1 LISA webpage: http://lisa.nasa.gov/ 


2 JWST webpage: http://www. jwst.nasa.gov/ 
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range, but in fact becomes increasingly less relevant for higher 
masses. 

It is also clear that deuterium shell burning is not going to 
occur close to the outer surface. While the characteristic tem- 
perature for a given mass shell scales as T « GM/R, we note 
that, from Eq. (10), GM/R œ t = M/M at times t > M/M. 
In this regime, the central temperature is almost spatially con- 
stant, consistent with the results of Hosokawa et al. (2012a), but 
already considerably larger than the deuterium burning temper- 
ature of ~10° K. Deuterium burning will thus occur when the 
first term in Eq. (10) is still relevant. From Fig. 1, it is evident 
that the radius of the shells changes considerably in that regime, 
while the mass in the shells is almost unchanged. Adopting a 
scaling relation of T « GM/R with almost constant M implies 
that the radius has to change by 2-3 orders of magnitude for 
the atmospheric temperature of 5000 K to increase to a value 
above 10° K. We therefore expect that the radius of deuterium 
burning will correspond to a fixed fraction of the protostellar 
radius as long as the protostar maintains its bloated envelope. 
From a comparison with Hosokawa et al. (2012a) at a stellar 
mass of 1000 Mo, we obtain the normalization of this relation 
to be 


1/3 
M 
Rp = 10 Ro (=) R 


© 


(A.2) 


As long as the protostar remains on the Hayashi track, the deu- 
terium burning shell will not be able to move toward the at- 
mosphere, and, therefore, it does not have an impact on the 
evolution of the protostar. 


A.2. The pp-cycle 


We consider first the energy production in the nuclear core via 
the pp-cycle, as the composition of the star is initially primordial. 
As shown by Hosokawa et al. (2012a), the typical temperature 
in the core is ~10° K, and the nuclear density ~1 g cm™?. We 
adopt here the expression of Padmanabhan (2000) for the energy 
production rate, noting that 


_ 24x 10*px? 22 


&p = ergs! el, (A.3) 


3 
E 
where To denotes the temperature in units of 10° K and X ~ 0.75 
the mass fraction of hydrogen. While the pp-cycle considered 
here has only a moderate temperature dependence, we note that 
the CNO cycle scales as T”? around temperatures of 10° K, and 
the triple-a process as T^ around temperatures of 10° K. For the 
latter cases, we will thus need to take into account the higher 
temperatures within the center of the core, as they substantially 
contribute to the energy production via nuclear burning. Now, 
with Gan ~ 1 gcm™ and T9 = 0.1, we have ee = 43 erg g`! el, 
The luminosity provided by the pp-cycle is then given as 


Lop = gn Me = 21.502/2 Lo. 


710 (A.4) 


For comparison, the Eddington luminosity of the protostar is 
given as 


Lea = 3.8 x 10* Mt Lo = 2.7 x 107r tno Lo. (A.5) 


The ratio of these luminosities is given in Fig. A.1, showing that 
it remains considerably smaller than unity for stellar masses up 
to 108 Mo. Equating the two expressions, we find that the lumi- 
nosity resulting from the pp-burning is relevant only at very late 
times tp) = 1.1 x 10% yr. 
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Fig. A.1. Ratio of nuclear luminosity provided via the pp-cycle 
(Eq. (A.4)) over the Eddington luminosity of the star (Eq. (A.5)) as a 
function of protostellar mass for different mass accretion rates. Higher 
accretion rates imply that less time was available for contraction, thus 
reducing the nuclear luminosities. 


For practical accretion rates of — 1073—10 and protostellar 
masses between 10°—108 Mo, the contribution from the pp-chain 
will never be relevant, and can thus be neglected for the overall 
evolution of the stars. 

We now estimate the amount of helium produced in the core. 
During one fusion event, an average energy of 4.3 x 1075 erg is 
released, implying a helium production rate of 


Been = zc (A.6) 
4.3 x 10> erg 

The helium mass in the core thus evolves as 

Mu, sn = 4mp X Nuie,pp- (A.7) 

An integration yields 

Muepp = 6.0 x 1078 51? Mo. (A.8) 


Equating the resulting mass with the total mass in the core, a 
helium core can be expected after a time of 1.2 x 10! yr. In 
summary, when only the pp-cycle is considered, nuclear burning 
is rather inefficient, and never going to have a significant impact 
on the evolution of the star. 
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